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Abstract 

Using the tomographic probability representation of qudit states and the inverse spin-portrait 
method, we suggest a bijective map of the qudit density operator onto a single probability distri- 
bution. Within the framework of the approach proposed, any quantum spin-j state is associated with 
the (2j + l)(4j + l)-dimcnsional probability vector whose components are labeled by spin projections 
and points on the sphere S 2 . Such a vector has a clear physical meaning and can be relatively easily 
measured. Quantum states form a convex subset of the 2j(4j + 3) simplex, with the boundary being 
illustrated for qubits (j = 1/2) and qutrits (j = 1). A relation to the (2j + l) 2 - and (2j + l)(2j + 2)- 
dimensional probability vectors is established in terms of spin-s portraits. We also address an auxiliary 
problem of the optimum reconstruction of qudit states, where the optimality implies a minimum rela- 
tive error of the density matrix due to the errors in measured probabilities. 

Keywords: spin tomography, spin portrait, qubit, qudit, probability representation. 



1 Introduction 

In the early years of quantum mechanics, it was proposed by Landau |1| and von Neumann [2] to 
represent the quantum states by the density matrices. This approach turned out to be applicable to 
spin states as well j.'jj . The density matrix formalism proved to be very useful and all physical laws 
such as the time evolution and energy spectrum were formulated in terms of this notion. However, 
many alternative ways to describe a spin-j state were proposed; for instance, with the help of different 
discrete Wigner functions (these and other analytical representations were reviewed in |4j) and a fair 
probability-distribution function w(m, n) called spin tomogram [5l|6]. The latter function depends on the 
spin projection m along all possible unit vectors n G S 2 . All these proposals are, in fact, merely different 
mappings of the density operator. Inverse mappings are also developed thoroughly and are based on the 
fact that, if a (quasi) probability distribution is given, the density matrix can be uniquely determined. On 
the other hand, there is a redundancy of information contained in spin tomogram w(m, n). An attempt 
to avoid such a redundancy was made in [7-15 . According to [7], the density matrix can be determined 



by measuring probabilities to obtain spin projection m = — j, — j + 1, . . . , j if a Stern-Gerlach apparatus is 
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oriented along Aj + 1 specifically chosen directions in space. In [12] it was shown that the density matrix 
can also be reconstructed if the probabilities to get the highest spin projection m = j are known for 
(2j + l) 2 appropriately chosen directions in space. In other words, one deals with the values of function 
w(j, rife), k = 1,2, ... , (2j + l) 2 and solves a system of linear equations to express the density matrix 

elements in terms of the probabilities w(j, n^). The conditions on vectors and an inverse 

method were also presented |12| . 

In this paper, we address the problem of identification of a qudit-j state with a single probability 
distribution vector. Moreover, such a vector must have a clear physical interpretation. These arguments 
make this problem interesting from both theoretical and practical points of view. An attempt to construct 
a bijective map of the density operator onto a probability vector with some interpretation of vector 
components was made in the series of papers 16 -20] (see also the recent review [2l]). The problem was 



shown to have an explicit solution for the lower values of spin j. 

In [22], a unitary spin tomography was suggested for describing spin states by probability distribution 
functions w(m, u), where u is a unitary (2j + l) x (2j + l) matrix. Information contained in this probability 
distribution function (called unitary spin tomogram) is even more redundant than that contained in spin 
tomogram w(m, n). Nevertheless, these redundancies can be used to solve explicitly the problem under 
consideration, i.e., to find an invertible map of the spin density operator onto a probability vector with 
a clear physical interpretation for an arbitrary spin j. The aim of our work is to present a construction 
of the following invertible map. It provides the possibility to identify any spin state with the probability 
vector "P with components V(m,Uk)- Here, random variables m and u k are the spin projection and 
unitary matrix, respectively. The spin projection takes values m = —j, — j + 1, . . . , j and a finite set 
of unitary matrices {uk}^=i describes the unitary rotation operations in the finite-dimensional Hilbert 
space of spin states. 

Thus, the function V(rn,Uk) is the joint probability distribution of two random variables m and Uk- 
The function V(m,Uk) is related to the spin unitary tomogram u>(m, u&) by the formula 

V(m,u k ) 

w{m,u k ) - 



This fact makes it possible to determine the density operator p. On the other hand, the probability 
distribution V(m, Uf~) contains extra information in comparison with that contained in the density matrix. 
The matter is that, if the density matrix p is given, a formula for the probability distribution V(rn, u k ) 
in terms of p does not exist. To obtain such a formula, one needs to take into account some additional 
assumptions. For example, we can assume the uniformity of the probability distribution of unitary 
rotations u k , i.e., each matrix u k (or direction in the case u G SU(2)) is taken with the same 
probability 1/N U , where N u is an appropriate number of unitary rotations. If this is the case, the density 
matrix p provides an explicit formula of the probability distribution V e q(jn, One can choose another 
nonuniform distribution of unitary rotations, but knowledge of this distribution is necessary information 
to be added to information contained in the density matrix. 

In the context of probability theory, the vector "P (or the point on a simplex determined by this 
vector) is defined by the set of spin projections m and unitary rotations Uk, which can be chosen with 
some probability. In view of this, the probability V(rn,Uk) under consideration is a fair joint probability 
distribution function of two discrete random variables. 

This paper is organized as follows. 
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In Sec. [2j a short review of spin and unitary spin tomograms is presented and a spin-s portrait 
method is introduced. In Sec. |3j we review a density matrix reconstruction procedure proposed in 12 



In Sec. [4j a representation of quantum states by probability vectors of special form (2l] is given. In 
Sec. [5j an inverse spin-portrait method is presented. This method allows to construct a map of the 
spin density operator p onto the probability vector "P. Two cases of used unitary rotations are given: 
u G SU(2) and u S SU(N), N = 2j + 1. In Sec. [6j the inverse map "P — > p is presented in explicit 
form for SU(2) rotations. In Sec. [7J particular properties of symbols V(m, n) are analyzed: star product 
kernel is presented and relation to symbols w(m, n) is considered. In Sec. [8j examples of qubits (j = 1/2) 
and qutrits (j = 1) are presented. In Sec. [9l the probability vector V is considered on the corresponding 



simplex and a boundary of quantum states is presented for qubits and qutrits. In Sec. 10, conclusions 
are presented. 



2 Unitary Spin Tomography and Spin-s Portrait of Tomograms 

We begin with some notation. Unless stated otherwise, qudit states with spin j are considered. Any 
state vector of such a system is uniquely determined through the basis vectors \jm), which are eigenvectors 
of both angular momentum operator J z and square of total angular momentum, i.e., J 2 = J 2 + J 2 + Jf . 
The spin projection m takes the values — j, — j + 1, . . . , j. 

A unitary spin tomogram of the state given by the density operator p is defined as follows: 

w^\m,u) = {jm\v) pu\jm) = Tr^p u\jm) (jm\vf\ = Tr^p U^\m,u)\ (1) 

where, in general, u is a unitary transform of the group SU(N). For the sake of convenience, starting 
from now we will identify u and its matrix representation u in the basis of states \jm), assuming that 
the matrix u defines, in fact, the unitary transform u. The tomogram w^\m,u) is a function of the 
discrete variable m and the continuous variable u. The operator U^> (m,u) = u\jm)(jm\v) is called the 
dequantizer operator because it maps an arbitrary density operator p onto the real probability distribution 
function w^\m, n). The dequantizer satisfies a sum rule of the form s }^m=~j U^\m, u) = I for all u. In 
view of this fact, the tomogram w^(m,u) is normalized, i.e., ^ 3 m= _jW^(m,u) = 1. 

The particular case u G SU{2) leads to the so-called spin tomogram w^'(m, n), where the di- 
rection n = n(8,4>) = (cos 4>s'm9, sin c^sin 6, cos 6) determines the dequantizer operator (m,n) 



R(n)\jm){jm\B) (n) (some properties of spin tomogram w^\m, n) were discussed in 23-25]). Here, we 
introduced a rotation operator R(n) defined through 

R(n) = e~ i(n ^ n ± = (-sine/., cos 0,0). (2) 

The inverse mapping of spin tomogram w^\m,n) onto the density operator p is relatively easily 
expressed through the quantizer operator D^\m,n) as follows: 

j 27r n 

p= ^f^f sindde w (j) {m,xi{6,cp))b {j) {m,xi{6,(t))). (3) 
m=-j o 

Different explicit formulas of both dequantizer and quantizer operators are known (5}[6}[26]-[28) , with the 
ambiguity being allowed for the quantizer operator. In this paper, we preferably use the orthogonal 
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expansion of the form 29 



2j 2j 



L=0 



L=0 

P(m,n) = ^(2L + l)^ ) (m)«(n)4 j) J Rt( n ) = ^( 2 L + l)/^(m)5?(n), 



(4) 



(5) 



L=0 



L=0 



where the coefficient f^\m) is an L-degree polynomial of the discrete variable m, and the operator 

S^\n) is the same polynomial of the operator variable R(n)j z R)(n) = (J • n). 
For instance, in the case of qubits (j = 1/2) and qutrits (j = 1), we have 

Zo (1/2 V) = ^, fi 1/2 \m) = V2m, S^\n) = ±i, s[^(n)- 



^ 1] (n) 



1 /, S[ l \n) = -^(J • n), S«(n) = i= fs(J • n) 2 - 2/ ) . 



V3 



v / 2 V 



1 

7q 



2(J-n), (6) 

(7) 
(8) 



In the case of an arbitrary spin j, /^(m) = l/\/2j + 1, fi\m) = ' yj j(j + l)(2j + 1), and all 

other coefficients (m) are expressed via the recurrence relation that was presented in 



fjj'(m), f^\(rn), and f^ 2 ( m )- Expansions (4) and (5) are orthogonal in the following sense: 



rC?) 



29 



and relates 



iv(^' ) (n)5g ) (n))= £ /^(m)/^(m) = ^, 



V), 



(9) 



m=-j 



This means that functions f^\m) are orthogonal polynomials of a discrete variable, with the weight 
function being identically equal to unity. Using the theory of classical orthogonal polynomials of a 

, it is not hard to prove that the function (m) is expressed through the 



discrete variable 



30 



31 



discrete Chebyshev polynomial t n (x,N) or Hahn polynomial h^'^\x, N) as follows: 
1 



/l V) = + m, 2j + 1) = ^4°' 0) (i + m, 2j + 1) 



d L 



' (2j + £ + !)! 
(2L + l)(2j-L)!' 



(10) 



2.1 Spin-s Portrait 

The tomogram w^(m, u) of a system with spin j is a function of the discrete spin projection m. This 
means that the tomogram w^(m, u) can be represented in the form of the following (2j + l)-dimensional 
probability vector: 

/ w^(j,u) \ 
w^\j — 1, u) 



V 



in 



(11) 



-Ji u ) j 
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We will refer to such a (2j + 1) vector as the spin-j portrait (in analogy with the qubit-portrait 
concept [32]). Note, that vector ( |11[ ) is a fair probability distribution vector since w(m,u) > and 
Em=— j w ( m i u) = 1 for all unitary matrices u. Since the components of vector are fair probabilities, 
the sum Eme.4 w^\m, u), where A C {m} J _-, can also be treated as a probability and has a clear physical 



(7) 

meaning. Summing some components of vector (11), one can construct a probability vector w s of less 



dimension (2s + 1), where s plays the role of pseudospin and can take values s = 1/2, 1, . . . , j — 1/2, j. 

(i) 

To be precise, w s reads 



w 



U) 



(«) 



( E 
E 



W 



IV 



m,u) 



m,u) 



me_4 2s +l 



(12) 



where *4fc C {m} 3 - for all 



1. 



,2s+l, U^f+Ufc 



{m}i-, and Anii = for all k^l. Vector (12) 



is referred to as the spin-s portrait of a qudit-j state. In the case s 
portrait of the form 



1/2, we obtain the so-called qubit 



IV 



IV 



0)| 

(i). 



m, u) 



m, u) 




meAl w^(m,u) 



1 -Y, m .eA2 wij) ( m > u ) 



J2rneA 2 w{3)( - m > U ) 



(13) 



with AiUA 2 = {m}{_ ,, ^li n ^l 2 



In the particular case of ^.1 = {j} and ^2 



portrait (13) reads 



w 



(j) 
1/2 



E 



0"). 



m, n 



1 — ioCj) (j 5 u) 



{m}-^ 1 , the qubit 



(14) 



Qubit portraits of this kind are implicitly used in several reconstruction procedures considered in subse- 
quent sections. The qubit-portrait method is introduced in [32] and successfully applied not only to spin 
systems [33] but also to the light states [34] . 

The inverse spin-portrait method is to construct a single probability distribution vector with the help 
of several spin-s portraits. Indeed, one can stack N s spin-s portraits into the following single 

probability vector of final dimension N s (2s + 1): 



1 



w^(m) \ 
(j) 



w 



(«a) 
U)( u Na ) J 



(15) 



This idea is elaborated for the case s = j in Sec. [5] Here, we restrict ourselves to the discussion of a 
particular case s = 1/2 (qubit portrait). Only one component of two- vector (13) contains information 



on the system. The density operator p is determined by (2j + 1) real numbers without regard to the 



normalization condition. This fact shows that, if a single probability distribution (15) contains complete 



information on the system, then it should comprise at least N±/ 2 = (2j + l) 2 different qubit portraits. 
Then the dimension of vector (15) is 2(2 j + l) 2 , but only (2j + l) 2 components are relevant. In the 
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following section, we review the reconstruction procedure that was suggested in |12| and that employed 
(2j + l) 2 qubit portraits of the form (14) with Uk G SU(2), k = l,...,(2j + l) 2 . The symmetric 
informationally complete POVM (positive operator-valued measure) to be outlined in Sec. [4] can also 
be treated as an implicit use of (2j + l) 2 qubit portraits (14) with unitary matrices Uk £ SU(N), 



N = 2j + 1, k = 
and (jj\u{ui\jj) 



1, . . . , (2j + l) 2 which satisfy additional requirements ^1=1^ ^k\ji){jj\u\ 
= l/(2j + 2) for all k ^ I. 



(2j + 1)/ 



3 Amiet— Weigert Reconstruction of the Density Matrix 



In this section, we review the approach [12] to reconstruct the density matrix of a spin j through 
the Stern-Gerlach measurements, where one measures the probabilities to obtain the maximum spin 
projection m = j for (2j + l) 2 appropriately chosen directions in space. 

According to ([l]), for each fixed direction n^, k = 1, . . . , (2j + l) 2 , the probability to obtain the spin 
projection m = j is given by the formula 



w®(j,n k ) = (jj\R\n k )pR(^k)\jj). 
We denote by W a vector comprising all these probabilities 

W - ( wOO (j, m) u>W (j, n 2 ) ... tub") (j, n {2 j+iy] 



(16) 



(17) 



Moreover, the (2j + l)x(2j + 1) density matrix p can be written in the form of a (2j + l) 2 -component 
vector 



Pn P21 



P2j+l,l Pl2 P22 



P2j+1,2 



Pl,2j+1 P2,2j+1 



P2j+l,2j+l ) , 



(18) 



which implies that the density-matrix columns are merely stacked up in a single column. Let U^(m, n) 
be a vector constructed from the operator U^' (m, n) using the same rule. Then the probability (m, n) 
is nothing else but the scalar product of two (2j + l) 2 -component vectors 



w (i) (m,n) = (\J ij) (m,n) ■ pj = Tr U (j) (m,n) T p 
From this, it is not hard to see that the map p —> W reads 

W = \\M\\p, 

where ||M|| is a (2j + l) 2 x (2j + l) 2 matrix of the form 

U^(i,n 2 ) T 



(19) 



(20) 



IMI 



V U«(j,n 



(2j+l) 2 



(21) 



/ 
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Reconstruction of spin-j density matrices: j = 1 (a, b) and j = 3 (c, d). The procedure 
is represented by a and c; it is based on fixing (2j + l) 2 specifically arranged direc- 



12 



Figure 1: 
developed in 

tions (23) and measuring the probabilities of the spin projection onto each direction to take the value j. 
The method used in [7] is illustrated by b and d; in this case, one needs to measure 4j + 1 spin-j portraits 
w^(rife), where are specifically chosen to form a cone. 



Whenever det ||M|| ^ 0, there exists an inverse map W — > p and 

p=\\M~ 1 \\W. (22) 

In ' 



12 



, it is shown that a particular choice of directions n^, k = 1, . . . , (2j + l) 2 ensures that the 
condition det ||M|| ^ is fulfilled. Moreover, such a choice of directions simplifies significantly the 
inversion procedure because it relies on the Fourier transform. Namely, the following directions were 
proposed: 

rifc = n qr = (sin 9 q cos (p qr , sin 6 q sin (p qr , cos 6 q ), < q, r < 2j, (23) 
with < 9 q < ir, 6 q ^ 9 q i if q ^ q' , and 

(p gr = 2n (r + qA), 0<A< — - — . (24) 

Examples of such a choice of the directions for spins j = 1 and j ' = 3 are shown in Fig. [TJ We see 
that the directions are divided into groups that form nested cones. Some modifications to free cones and 



spirals were presented in 15 . An arbitrary choice of the directions was discussed in [14| 



4 Quantum States as Probability Distributions 

To begin, we showed in Sec. [2] that any quantum state can be interpreted as a probability distribution 
function. Since tomogram ([!]) is a function depending on the unitary matrix u (direction n in the case of 
u € SU(2)), there is a redundancy of information contained in the tomogram. It is tempting to reduce 
such a redundancy and associate a quantum state with the single probability distribution. The proposal to 
associate quantum states with single probability vectors was made in |16f[21| . To avoid any redundancy, it 
was suggested to use the minimum informationally complete POVM (positive operator- valued measure). 
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This suggestion is related to constructing the minimum tomographic set discussed in [36] . Using the 
language of spin states, any spin-j state is associated with the following (2j + l) 2 -component probability 
vector: 

/ tr(p^) \ 

VP ' , (25) 
\ tr(p-E (2 j + i)2) J 

where Ei, i = 1, . . . , (2j + l) 2 are the corresponding POVM effects. There are many ways to introduce the 
minimum informationally complete POVM effects, and a possible choice that is valid in any dimension 



was presented in 19 . A relatively new tendency is to use the symmetric informationally complete (SIC) 

POVMs. If this is the case, the effects {Ei}^ 1 ^ 
condition: 

' if k^l. (26) 



As stated in 



20 



, the effects {Ei\\ 



are one- dimensional projectors satisfying the following 
1 



(2j + l) 2 (2j + 2) 



1 



that meet the above requirements were found numerically in 
dimensions (2j + 1) < 67 and analytically in dimensions (2j + 1) = 2 — 15, 19, and 24. 

It is worth noting that the Amiet-Weigert construction considered in the previous section can be 



treated as a single probability distribution. For this, one needs to normalize vector (17). In this case, 



the mappings (20) and (22) are slightly modified 



w 



I Ml 



IMIHW 



v (2j+i) 2 



\M\\ ik pk 



^2j+l v-(2j+l) 2 
Zui=l Z^fc=i 



I^H(i-l)(2j+l)+i,fe^ / fc 



(27) 



5 Inverse Spin-Portrait Method 

In this section, starting from spin tomograms, we associate each quantum state with the corresponding 
probability distribution vector, which has a clear physical meaning and can be measured experimentally. 

The unitary spin tomogram tow) (jn, u) is a function of the unitary matrix u and depends on a discrete 
parameter m = —j, . . . , j— 1, j. Fixing the unitary rotation Uk, we obtain a (2j + l)-component probability 
vector Wj (uk), called the spin-j portrait of the system (see Sec. 2.1). Given only one spin-j portrait 
of the system, it is impossible, in general, to define without doubt a state of the system. Nevertheless, 
the state is determined if one has an adequate number N u of different spin-j portraits. Then we can 
introduce a joint probability distribution function of two random variables m and u k 



V(m,Uk), m = -j, -j + 1, . . . , j, k=l,2,...,N u . 



(28) 



The physical meaning of this joint probability distribution is that, if one randomly chooses a unitary 
rotation from the set {uu\k=i an( ^ a s P^ n projection m within the interval —j<m<j, the value of 
V(m, itfe) gives the probability of the detector's click. Function (28) can also be written in the form of 



the following N u (2j + l)-component probability distribution vector: 



r(j, Ul ) 



V(-j,ui) V(j,u 2 ) 



r(-j,u 2 ) 



'Pi- 3, u Nu ) 



(29) 
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with the normalization condition of the form 



m,u k ) 



(30) 



-3 k=l 



Since the tomogram w^\m,Uk) 1S nothing else but the probability to obtain the spin projection m 
if the rotation Uk is fixed, the relation between the spin tomogram and the joint probability distribution 
function V{m^Uk) reads 



w^'(m, Uk) 



i— dm 



(31) 



where the denominator has the sense of the probability pk to choose the unitary rotation u^. If the 
probabilities {pk\k=l are known a priori, then the vector "P is easily expressed via spin-j portraits 
namely, 

/ Piwf(m) ^ 



(32) 



If the unitary rotations «&,&=!,..., are equiprobable, then 



eq 



1 

A^ 



wr'(ni) 



w 



(«2) 



(33) 



It is worth mentioning that, even if a priori the probabilities {pk}^=i are n °t known, formula (31) 



provides a direct way of mapping V onto the vector "P e q- 

Let us now consider an open problem of the minimum number N u of spin portraits. In other words, 
N u is the number of unitary rotations Uk that is needed to identify any quantum state with a single 



probability vector of the form (32) and to minimize the redundancy of information contained in this 



vector. Subsequently, two main cases are presented, namely, the use of SU{2) rotations and SU{N) 
rotations with N = 2j + 1. These particular problems can be of interest for experimentalists because 
SU(2) rotations can be relatively easily realized in some modifications of the Stern-Gerlach experiment, 
while SU(N) rotations may require more difficult apparatus. On the other hand, it will be shown that, 
to extract information on the system, one can use a smaller number of SU (N) rotations than in the case 
of SU (2) matrices. 



5.1 SU(2) Rotations 



Like the Amiet-Weigert scanning procedure ( 20 ) , the map of the density operator p onto the probability 

V = IIQUp, (34) 



vector (32) can be written as follows: 



9 



where ||Q|| is an N u (2j + 1) x (2j + l) 2 rectangular matrix of the form 



IIQII 



I Pl uO)(j, ni ) T \ 
Pl uW(j-i, ni ) T 



PNu \j( j) U,n Nu ) T 
p Nu V^(j -l,n Nu ) T 



\ PN, 



U(i)( 



'3i njv u 



(35) 



The map ( |34[ ) is invertible iff rank||Q|| = (2j + l) 2 . Since the rank of a matrix is equal to the number 
of linearly independent rows, the number N u can be defined as the minimum natural number such that 
the set {"U"( J )(m, n^)}, m = k = 1,...,N U contains (2j + l) 2 linearly independent vectors. 

, each vector U^(m, n) can be resolved to the sum of orthogonal vectors S^(n), 



29 



According to 
L = 0, . . . , 2j; namely, 



2.7 



uW(m,n) = ^/l j) (m)S^(n), 



id), 



(36) 



L=0 



where ( S?(n) • S#(n)J = Em=-j /^M/i'V) = $ll 



r(j), 



erator S'^(n) acting on the Hilbert space of spin-j states (see Sec. 2). The operator 5^''(n) is shown 



The vector (n) corresponds to the op- 



to be the same polynomial f^\m) of degree L, with the argument being replaced: 
Jxn x + Jy n y + Jz n z- Suppose L = 0, then there exists only one linear independent vector of the form 
^ii=o( n ) wmcn corresponds to the identity operator /. If L = 1, no more than three linear independent 
vectors S^4i k = 1, 2, 3, which correspond to the operators (J-ni), (J-i^), and (J -113), respectively, 

can exist. Note that the three vectors S^2=i ( n A;) 5 k = 1,2,3 are independent iff vectors £ M 3 are not 
coplanar, i.e., their triple product (ni • [112 X n3]) 7^ 0. Taking into account the normalization conditions 
n| = 1 and J 2 = j(j + 1)7, in the case L = 2, we obtain five linear independent vectors S^=2( n k)- Using 
the matrix representation of the operator S^\n), we see that it is composed of 2L + 1 independent 
/-diagonal operators (1 = for a diagonal one, I = 1 for a super-diagonal one, / = — 1 for a sub-diagonal 
one, and so on for all —L < I < L). Increasing L by unity, two more diagonals are filled. We draw the 
conclusion that for a fixed L the maximum number of linearly independent vectors s9 (n*.) is equal to 



m 



(J-n) 



2L + 1. Since vectors (n^) and S^y (n^) with different L and L' are orthogonal, the total number of 
linear independent rows Uj(m, n) equals 1 + 3 + 5 + • • • + N u = (N u + l) 2 /4. On the other hand, it must 
be equal to rank 1 1 Q 1 1 = (2j + l) 2 . From this, it is readily seen that N u = Aj + 1. 

The directions n^, k = 1, . . . , 4j + 1 cannot be chosen arbitrarily because of the condition rank||Q|| = 
(2j + l) 2 . As was shown above, the directions {nfc}^^ 1 are divided into sets of one, three, five, and so 
on directions. Without loss of generality, it can be assumed that these sets are {ni}, {n^}| =1 , {nfc}| =1 , 
. . . , {n/j}^^ 1 , respectively. If this is the case, the requirement rank 1 1 Q 1 1 = (2j + l) 2 is equivalent to the 
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condition 



where A q , q = 1, 



A X A 2 • . . . • A 2j ^ 0, (37) 
,2j are expressed through = (sin0fc cos ipk, sin 9k sin cp^, cos 6*^) and associated 



Legendre polynomials P, (a;) as follows: 



/ P g (O) (cos0i) 



A„ = det 



P 



(o), 



COS #2 J 



P 9 (O) (cos0 3 ) 



P g (O) (cos0. 



2qj 



In the particular case of q = 



Pq 9 ^ (cos 9\ ) cos qip\ 

Pq q \cOS 9 2 ) COS<?V?2 
Pg 9) (cOS 3 ) COS ^3 

P^(cos 9 2q ) cos qip 2q 

Ml 



P^ (cos 0i ) sin gt^i 
Pj^(cos 2 ) sing^ 
P^cos 3 ) sing^ 3 



P g (9) (cos 2(? ) sin qip 2q 
Mi 



P g W; (cOS0 2 g+l) COS^ 2 g+l Pg (COS 2(?+ l ) Sm g^ 2 g+l / 

1, we have 



/ P 1 (O) (cos0i) P 1 (1) (cos0i)cos9?i P 1 vi; (cos0i)sin<^i 
P 1 (O) (cos0 2 ) P 1 (1) (cos0 2 )cos( / 9 2 P 1 (1) (cos0 2 )sin( / 92 
V P 1 (O) (cos0 3 ) P 1 (1) (cos0 3 )cos9? 3 P 1 (1) (cos0 3 )sin 



Ai = det 



(ni • [n 2 x n 3 ]). 



(39) 



<P3 



It is worth mentioning that there exists an optimum choice of the directions {nk}k=i that provides 
minimum possible errors of the reconstruction procedure due to errors in measured probabilities V(m, n&). 
Indeed, according to the results of computational mathematics (see, e.g., [35]), the errors of the vector p 
defined by formula (34) are directly proportional to the condition number fi(Q) of the matrix ||Q||. The 
greater the product Ai • . . . ■ A 2 j, the smaller fi(Q) and, consequently, the errors of the reconstruction 



procedure. For qubits (j 



1/2), the optimum choice is three orthogonal vectors {nfc}| = 



because their 



triple product takes the maximum value in this case. As far as higher spins are concerned, maximization 



of expression (37) is performed numerically, and the optimum directions {n^jj^^ are shown in Fig 



A particular case of 6% = 9 2 = . . . = 0<y+i = corresponds to the Newton- Young reconstruction 
procedure |7j and implies the limitation P^ m \cos9) ^ for all < m < L, 1 < L < 2j. A schematic 
illustration of this reconstruction procedure is given in Fig. [TJ 

An alternative way to meet the requirement rank||Q|| = (2j + l) 2 is to ensure linear independence 
of vectors {S^(rifc)}^^ 1 for all L = l,...,2j. Linear independence of these vectors is equivalent to 



nonzero Gram determinant det 
we obtain 



Mi 



det 



det 



'Lk=l 



P 2 \ n i ■ n k, 



2L+1 



,fc=l 



i,fc=l 



Since I8$>(m) ■ (n fc )) = P^(n t • n fc 



Ml 



det 



if (Pi-**) 



4j+l 
i,k=l 



^0. 



(40) 



In the case of qubits, one has 



Ai = det 



< (ni • ni) (ni • n 2 ) 
(n 2 • ni) (n 2 • n 2 ) 
\ (n 3 • ni) (n 3 • n 2 ) 




(m • [n 2 x n 3 ]) V 0. 



(41) 
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Figure 2: Optimum choice of the directions {n^},^ that ensures the minimum possible errors of the 
density operator p with respect to the errors of the probability vector "P, where j = 1/2 (a), j = 1 (b), 
j = 3/2 (c), j = 2 (d), j = 5/2 (e), and j = 3 (f). 



As far as qutrits (j = 1) are concerned, we obtain 



A' r det 



, 3(n r n 2 ) 2 -l 3(n r n 3 ) 2 -l 3(n r n 4 ) 2 -l 3( ni -n 5 ) 2 -l \ 

2 2 2 2 

3(n 2 -ni) 2 -l -. 3(n 2 -n 3 ) 2 -l 3(n 2 -n 4 ) 2 ~l 3(n 2 -n 5 ) 2 -l 

2 1 2 2 2 



3(n 3 -m) 2 -l 3(n 3 -n 2 ) 2 -l 



1 



3(n 3 -n 4 ) 2 -l 3(n 3 -n 5 ) 2 -l 



2 2 

3(n 4 -ni) 2 -l 3(n 4 -n 2 ) 2 -l 3(n 4 -n 3 ) 2 -l , 

2 2 2 L 

3(n 5 -m) 2 -l 3(n 5 -n 2 ) 2 -l 3(n 5 -n 3 ) 2 -l 3(n 5 -n 4 ) 2 -l 

2 2 2 2 



3(n 4 n 5 ) 2 -l 
2 

1 



^0. 



(42) 



/ 



5.2 SU(N) Rotations in Hilbert Space 

If the matrix u in the definition of spin tomogram w^\m,u) is an element of the group SU(N) with 
2 < N < 2j + 1, the tomogram w^\m,u) is also referred to as a unitary spin tomogram 22 . The case 
N = 2j + l corresponds to the most general form of unitary rotations in the Hilbert space of spin j. By the 
previous statement, we can assume linear independence of vectors {U(m, Uk)}, where the spin projection 
m takes values m = —j, — j + 1, . . . , j — 1 and k runs from 1 to N u . Here N u denotes the minimum number 
of spin portraits w^(ufc), G SU(N), N = 2j + 1, that are needed to construct a bijective map. Since 

I for all k = 1, . . . ,N U , we have one more independent vector. Thus, the total 



3 

m,u k ) 
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number of linear independent vectors {U(m, u^)} equals 2jN u + 1. On the other hand, this number must 
be equal to (2j + l) 2 for the map to be invertible. Consequently, N u = ((2j + l) 2 — 1 ) / 2 j = 2j + 2. 

We see that one needs fewer SU(N) rotations than SU(2) rotations in order to map a quantum state 
onto a single probability distribution. Nevertheless, this advantage is accompanied by the complexity of 
the experimental realization of SU(N) rotations with regard to 577(2) rotations in the Hilbert space. For 
intermediate cases 2 < N < 2j + 1, we have 2j + 2 < N u < 4j + 1. 

Let us consider the case of SU(N) rotations with N = 2j + 1 in detail. As in the previous section, 
any quantum state is mapped onto a single (2j + l)(2j + 2)-probability vector "P as follows: 



V= \\R\\p, 

where the (2j + l)(2j + 2) x (2j + l) 2 rectangular matrix \\R\\ reads 



(43) 



/ 



\R\\ 



Pi 



\ 



PiUW(j-i,m) Q 



P2j+2U (j) (j, U 2 j+2) T 

P2i+2Utf)(j - l,n 2j+2 ) 

V P2i+2U (j) (-J,lt 2 j + 2)' r / 

Taking into account the condition rank||i?|| = (2j + l) 2 , we express the inverse map "P 
pseudo-inverse matrix \\R~* 



(44) 



37 as follows: 
p= \\R + \\1 



R\ 



-i 



\\tf\\V. 



p through a 



(45) 



The requirement rank||i?|| = (2j + l) 2 is equivalent to the linear independence of vectors U^'(m, itfe), 
m = —j, — j + 1, ... ,j — 1, fc = 1, 2, . . . , 2j + 2 and the vector I = Ylm=-j ^fc), which does not rely 

on /c. The linear independence of these vectors is achieved whenever the corresponding Gram determinant 
is nonzero. That yields the constraint of the form 



det 



/ 2 j + 1 


11 ... 1 


11 ... 1 




11 ... 1 


l 
i 

— i— 




A(uii4) 




A(ui4 i+2 ) 


i 
i 


A(u 2 u\) 


hj 




A{u 2 u\ j+2 ) 












i 

l T 


A(u 2 j + 2u\) 


A(u 2j+ 2u\) 







(46) 



where the Gram matrix is composed of blocks: I 2 j is the 2j x 2j unity matrix and A(uj.v/y) is a 2j x 2j 
matrix with elements ||A(t4jfcn^,)||j n = \ {jl\uku\,\jn)\ 2 . Using the orthogonal expansion of dequantizer (|5|) , 
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we rewrite the condition obtained in terms of vectors (itfc) as follows: 



T' = det 









A'(ui,U 2 j +2 ) 


A'(u 2 ,u x ) 


hi 




A'(u 2 ,U 2j+2 ) 










\ A'(u 2i+ 2,ni) 


A'(u 2j+2 ,u 2 ) 




hi 



(47) 



where A'{u k , uy) is a 2j x 2j matrix with elements ||A'(iifc, Uk')\\LL' 



iCi), 



Uk, 



V), 



. Note that V is 



nothing else but the volume of the 2j(2j+2)-dimensional parallelogram with edges S£ (u k ), L = 1, . . . , 2j, 
k = 1, . . . , 2j + 2. Moreover, < T' < 1 since all vectors S^(itfc) are normalized. 

The optimum choice of unitary matrices u k & SU(N), N = 2j + 1, follows the same line of reasoning 
as in the case of SU(2) rotations. Indeed, if the probability vector "P is measured experimentally within 
the accuracy <5"P, formula (45) yields the vector p defined with an error bar of bp. It is known that 

1 1 S^P 1 1 ™~ 

— ^ ||-p||-f ' w here p, is the condition number of the matrix (|47|) and || • \\ 2 is the Euclidean norm of a 



IIPII2 



l+v / l^T 



vector. It can be shown that p < -_ j_ y 
reconstructed state p. 



j. Consequently, the greater T' (or T), the less erroneous is the 



6 Inverse Mapping of a Probability Vector onto the Density Matrix 

In this section, we give an explicit expression of the density operator p of the system with spin j in 
terms of the single probability vector "P, which could itself be treated as the notion of quantum state. 
We consider distributions "P obtained by SU{2) rotations. It was shown in Sec. [5] that any vector "P is 
readily transformed into the vector "P e q- For this reason, we will focus attention on the map Peq —> P- 

Let us now recall that the direct map reads 

V eq (m,n k ) = ^±— w (i)( m ,n k ) = — -Tr(/5 £>^(m,n fc )), (48) 
where U^\m,n k ) is the dequantizer operator that can be resolved into sum This means that 

V cq (m,n k ) = --^(m.m) = — - £ /?(m)Tr (p§$\n h j) = ^Tr(p uf{m,n k )), 

(49) 

where we introduce the L-dequantizer operator W? (m, n k ) = (4j + f^' (m)Sj? (n^). 

Since the direct map is linear, it can be assumed that the inverse map is also linear, that is, 

43+1 j 

P= Yl 7 , eq(m,n fc )pW(m,A;), (50) 

k=l m=—j 

where P^) (to, k) is the quantizer operator to be determined. We already know that it is convenient to 
rearrange directions {n^}^ 1 and consider sets {n^}^^ 1 , L = 0, 1, . . . , 2j. In view of this fact, we treat 
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a solution of the form 



2j 2L+1 j 

^ = EE ^2 Pe q (m,n k )V^(m,k). 

L=0 k=l m=—j 



(51) 



If L-quantizers T>^\m, k) are known, we immediately have T>^\m, k) = ^ T>^\m, k). 



Proposition. The L-quantizer is expressed through operators S^\n k /) an d Gram matrix \\^{L)\\ 



L: (k~l)/2<L<2j 



whose matrix elements are \\^ \L)\\ kk i = Tr ^^(n^S^^n^ 
as follows: 

2L+1 

V ( i\m,k) = (Aj + l)f^(m) W^m^S^M, 



S?(n fc )-Sf (n fc /)) = P L (n k -n k . 



U) 



(52) 



k'=l 
2L+1 



with the operators {X^fcWL 1 \\^ 1 (-^ J )\\kk'S^ \n k i)\ forming a dual basis for the given basis 
I ) k=l 

□ 



Let us check that formula ( 52 ) gives an adequate solution of the problem. 

If the directions {n^.}^^ 1 are chosen properly and the requirement (37) is satisfied (which is equivalent 



to niii det \\JV(L)\\ £ 0), then any density operator p is resolved into a sum of orthogonal operators 

2j 2L'+1 

/ W 1 t \ Afl'l , 

in*/). 



(53) 



L'=0 k'=l 



Substituting formula (53) for p in (49), we obtain 



2j 2L'+1 

Tr(pS?(n fc )) = E E a(L',k')Tr[s^(n k/ )S { i\n k )) = £ a(L, fc')IW)ll 

L'=0 fc'=l 

^ jO') "\ \ IT 7 / X II SS/TW, 

\kk 



2L+1 



k'=l 



kk 1 



2L+1 



(54) 
(55) 



L=0 



fe'=l 



After combining (52) and (55), direct calculation of the right-hand side of Eq. (51) yields 

2j 2L+1 j 

E E E ^e q Kn fc )p[ j) (m,/c) 

L=0 k=l m=~j 

2j 2L+1 2L"+1 r j 

= E E E E 

L,L"=0k,k"=l k'=l l m=-j 



m) 



a (L", k')\\^ (L") \\ kk > ||^ _1 (L) \\ kk "S^ (nfc/) 



2j 2L+1 



■2L+1 



E E a ( L ' fc ') E ll-^WIUfc'll-^" 1 ^)!!^" ?(%) = ££a(L^')fK 



2j 2L+1 



L=0fc' fc"=l 



fc=l 



L=0 fc'=l 



5fc/i 



(56) 
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This concludes the proof. 
■ 

In fact, the proof above is followed by the relation between the L-dequantizer and the L'-quantizer 
TV iuf (m, n k )V { i) (m', n k ,)) = ff {m)jf (m')5 LL ,8 kkl . (57) 



To summarize the results of this section, we write the explicit form of the quantizer 



t>^(m,k)= 

L: {k-l)/2<L<2j 



SO')/ 



2L+1 



VY'(m,k) = (4J + 1) Yl fL( m ) E \\^-\L)\\kk>§r(nv). 



L: (k-l)/2<L<2j 



k'=l 



(58) 



7 Star-Product and Intertwining Kernels 



Suppose Veq i(rrii,n k J is the symbol (48) of a state pi, i = 1,2; then the operator pife is associated 



with a symbol V e q 3(^3,11^3) which is called the star product [38 - 40 of symbols V CC[ i(mx,n^ 1 ) and 
'Peq 2(^2, n-k 2 ) an d denoted by (V e q 1 * T^eq 2) {jnz-, nfc 3 )- Combining (49) and (58), it is not hard to see 
that 



4j+l 



T'eq 3K,nfc 3 ) = E E ^ (5 H m 3, n fc3 , m 2 , n^ 2 , mi, nfcjPeq 2(^21^2)^ 1 (mi, n fcl ), 

fci,fc2=l mi,ni2=-j 



where the star-product kernel K^'{m^, nfc 3 , m-2, nfc 2 , mi, n^) reads 

if ^ (m 3 , n fc3 , m 2 , n fe2 , mi , n fel ) = Tr ') (mi , n fcl )£> (j) (m 2 , n fca (m 3 , n fc3 

= (4i + l) £ E E /ff W/ff W/g(m 3 

{L i: (fc 1 -l)/2<L 1 <2j}{L 2 : (fc 2 -l)/2<L 2 <2j} L 3 =0 
2Li+12L 2 +l 



(59) 



(60) 



Let us recall that we have considered previously two maps of the density operator p onto the proba- 
bility distribution functions, namely, the map onto tomograms w^\m, n) depending on the continuous 
variable n = n{9,ip) G S 2 and the map onto the single probability distribution V cq (m,n k ) depending on 
the discrete variable n k , k = 1, 2, . . . , Aj + 1. Since both maps are invertible, symbols w^\m, n) and 
V eq (m,n k ) are related by intertwining kernels. Indeed, combining formulas ^ and (48), we obtain 



V cq (m,n k )= Y ^K^\ v (m,n k ,m\ii)w^\m,n) 

m'=-j S 2 



(61) 
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where K ~!_.<p(m, n^, m', n') = (4j + l) x Tr (i5W(m',n')t/'^(m, n^)). Using expansions (4) and (5) along 



with the orthogonality property Tr ( Sx/(n'),Sx(nfc) ) = 5lliPl{to! ■ n^), we arrive at 



(62) 



K { Jl v (m,n k ,m',n') = (4j + l)" 1 ^(2L + l)f { ^ (m')f^ (m)P L (n' • n fc ). 



2.7 



0')/ 



L=0 



In view of the same argument, the tomogram w^\m, n) is expressed through the joint probability 
distribution V eq (m' , n^/) as follows: 



4j+i j 

ij) (m,n)=^ ^ K^\ w (m,n,m, k')V cq (m' ,n k ,), 

k'=l m'=—j 



(63) 



where K^}^ w (m, n, m', if) = Tr (pW^fc')LK^( m, n) I . Taking into account explicit formula &58ti, we 
obtain 

2L+1 

K^Jm,n,m',k') = (4j + l) £ /jf'VO/rV) J] ll^WIU'^n* • n). (64) 

L: (k'-l)/2<L<2j k=l 



8 Examples: Qubits and Qutrits 

In this section, the results of the previous sections are specified for two particular cases of lowest spins, 
namely, qubits (j = 1/2) and qutrits (j = 1). As far as qubits are concerned, only the SU(2) rotations are 
possible. A quantum state is associated with the six-dimensional probability vector "P with components 
V(m, nfc), m = ±1/2 and k = 1,2,3. In other words, the probability vector "P is composed of three 
qubit portraits defined by directions {nfc}? =1 . If these directions are equiprobable (chosen with the same 



probability p k = 1/3), the corresponding probability vector is denoted as "P cq . Note that formula (31) 



defines the mapping "P — > T-'eq for any vector "P. For the map p — > P to be invertible, the limitation 
(ni • [ri2 x 113]) 7^ is imposed. The least erroneous reconstruction procedure (see Fig. [2]) takes place if 
all three directions are orthogonal, i.e., (ni • [112 x 113]) = ±1. In general, the inverse map (51) of the 
probability vector "P oq onto the density operator p reads 

1 3 

P = Peq( + l/2, ni) + Peq(-l/2, Hi)] I + £ ^eq( + l/2, Ilfc) - P C q(-l/2, life)] (<T • lfc), (65) 
^ fc=l 

where <x = (a^, a y , a z ) are the Pauli operators, and the vectors = Y^k'=i = l)IUfc' n fc'- Direct 

calculation provides 



li 



[n 2 • n 3 ] 



(ni • [n 2 x n 3 ])' 



h 



[n 3 • ni] 



(ni • [n 2 x n 3 ]) ' 



h 



ni • n 2 J 



(ni • [n 2 x n 3 ]) ' 



(66) 



Thus, the vectors {lfc}| =1 form a dual basis with respect to the directions {nfc/}|, =1 , i.e., (1^ • ny) = 8 k k'- 
This dual basis can be used to construct dual symbols of operators (41 42 . Note that vectors {lfe}| =1 are 
no longer normalized. Figure [3] illustrates the duality of basic sets {njfc}| =1 and {lfe}| =1 and, consequently, 
the duality of mappings p — > "P and "P — > p. 
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Figure 3: Duality of basic sets {n.k} k= i and {lfc}| =1 . Any qubit state is identified with a six-dimensional 



probability vector "P, which is composed of three spin-1/2 portraits w^Y ( n l)> Wy^' (112), and w ^2^(113) 
given by vectors ni, 112, and n3, respectively. Vectors {lfc}| =1 form a dual vector basis that determines 



the inverse map (65) of the vector "P onto a qubit density operator. The map is bijective whenever 



vectors ni, ri2, and n3 are noncoplanar. 



Using the special properties of Pauli matrices, one can easily calculate the star-product kernel (60). 
The result is 

K^ /2 \m 3 ,n k3 ,m 2 ,n k2 ,mi,n kl ) = 3| -S k3 ,xS k2tl + <W m 2«i-i(lfc 2 ■ n kl ) + 5 k2 ,im 3 mi(l k:i ■ n kl ) 

+m 3 m 2 (l k3 • lfca) + 2im 3 rri2mi([l k3 x \ k2 ] ■ n&J j. (67) 
The intertwining kernels (62) and (64) take the following form: 



K^l(m, n k , m, n') = i + 2m'm(ri ■ n k ), K^l(m, n, m , k') = 3 | ^k',i + 2m'm(lk> ■ n)| 



(68) 



As far as qutrits are concerned, any quantum state can be associated either with the fifteen-dimensional 
probability vector (34) parameterized by five SU(2) rotations or the twelve-dimensional probability vec- 
tor (43) written in terms of four SU(3) rotations in the Hilbert space. The former case implies that 
the vector "P comprises five qutrit portraits, each defined by the direction n^, h = 1, . . . , 5. The density 
operator is uniquely determined whenever these directions satisfy the condition (42). The latter case of 
SU(3) rotations implies the limitation (47) on unitary matrices u k i, k' = 1, . . . ,4. 
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The explicit formula of the density operator in terms of experimentally attainable probabilities 
V eq (m, ilk) reads 

1 1 3 

P = o [^eq(+l, ni) + P eq (0, 111) + 7>eq(-l, ni)] / + - [P C q(+l, n*) " Peq(-1, n*)] (J • lfc) 

fc=l 



+ - 



P cq (+l,n 2 ) 
P eq (+l,n 3 ) 
P eq (+l,n 4 ) 
\ 7%(+l,n 5 ) 



2P cq (0, ni ) + Peq(-l,ni) \ 
2P cq (0,n 2 ) + 7%(-l,n 2 ) 

2P e q(0,n3)+Peq(-l,n 3 ) 

2P eq (0,n 4 )+P eq (-l,n 4 ) 
2p eq (0 ! n 5 )+P eq (-l,n 5 ) ) 



1 



3(n r n 2 ) 2 -l 3( ni -n 3 ) 2 



3(n r n 4 ) 2 -l 3(n r n 5 ) 2 



2 

3(n 2 -ni) 2 -l , 
2 1 
3(n 3 -m) 2 -l 3(n 3 -n 2 ) 2 -l 
2 2 

3(n 4 -n 2 ) 2 -l 
2 

3(n 5 -n 2 ) 2 -l 



3(n 4 -ni) 2 -l 
2 

3(n 5 -ni) 2 -l 



3(n 2 -n 3 ) 2 -l 
2 

1 

3(n 4 -n 3 ) 2 -l 



3(n 2 -n 4 ) 2 -l 
2 

3(n 3 -n 4 ) 2 -l 
2 

1 

3(n 5 -n 3 ) 2 -l 3(n 5 -n 4 ) 2 -l 



-1 



3(n 2 -n 5 ) 2 -l 
2 

3(n 3 -n 5 ) 2 -l 
2 

3(n 4 -n 5 ) 2 -l 



-1 



/ 3(J 


rn) 2 


-27 \ 


3(J 


n 2 ) 2 


- 21 


3(J 


n 3 ) 2 


- 21 


3(J 


n 4 ) 2 


- 21 


V 3(J 


n 5 ) 2 


-21 ) 



(69) 



A similar reconstruction of qutrit states is proposed in [44 1. Higher spins (S = 4) are reconstructed 



in 



45 . The advantage of the proposed inverse mapping (plj) is that it can be applied to a system with 



an arbitrary spin j and provides the explicit solution in an operator form. 



9 Quantum States on 2j(4j+3)- Simplex 



We already know that any quantum state can be associated with the probability-distribution vector 
T> . If SU(2) rotations underlie the construction of the vector "P, this vector comprises (2j + l)(4j + 1) 
components V(m, n&). Consequently, it is represented by a point on the simplex with dimension (2j + 
1) (4j + 1) — 1 = 2j(4j + 3). Conversely, not all points on the 2j(Aj + 3) simplex can be associated with 
quantum states. Indeed, the condition p > is to be satisfied. Let us reformulate this requirement in 
terms of components 'P(m, n^). 



In view of the inverse map (51), we readily obtain the following condition: 



2j 2L+1 j 

EE E 

L=0 fc=l m=—j 



V eq (m,n k )f¥\m) £ W^ 1 (L)\\ kk , M > 0. 



2L+1 



(70) 



This implies that the matrix of the operator (70) in the basis of states \jm) is nonnegative. Nonnegativity 



of such a matrix is easily checked by Sylvester's criterion 43 ; namely, it is necessary and sufficient that 



all principal minors of matrix (70) are nonnegative. 

Let us consider the case of qubits (j = 1/2) in detail. 



19 




Figure 4: Convex subset of quantum states on the 2j (4j + 3)-simplex. In the case of qubits, 
P eq (+l/2, ni ) = pS 1/2) , 7> oq (-l/2,m) = l/3-PS 1/2) , P cq (+l/2,n 2 ) = P< 1/2) , P eq (-l/2,n 2 ) = 
1/3 - P(, 1/2) , P eq (+l/2,n 3 ) = P 3 1/2) , and P eq (-l/2,n 3 ) = 1/3 - P 3 1/2) (a). So any point on the V- 
simplex is uniquely determined by a point (^P^ 2 \p^ 2 \p^ 2 ^ inside a cube {0 < Pjj, 1 ^ < 1/3}| =1 , 
where quantum subsets are shown for (ni • [n 2 x n3]) = 1 (1), 0.44 (2), and 0.02 (3). In the 
case of qutrits, simplex is fourteen-dimensional (b), where we fix five directions in such a way that 
(ni • n 2 ) = (n 2 • n 3 ) = (n 3 • m) = (ni • n 4 ) = (n 2 • n 4 ) = (m • n 5 ) = (n 3 • n 5 ) = l/\/3. De- 
note Pj^ = V C(l (+l,iik), k = 1,2,3, then the cut set 'P cq (+l,n4) = V cq (+l,n^,) = 1/15, and 
V e q(— l,nfc/) = 1/15, k' = 1,...,5 corresponds to 1, the cut set 7^(4-1, iu) = 7 3 eq (+l, ns) = 1/15, and 
V e q(— 1, nfc') = 1/15, k' = 1, . . . , 5 corresponds to 2, and the cut set V eq (+l, 114) = "P C q(+l, ns) = 1/20, 
and 7 3 eq(— 1, n^/) = 1/20, k' = 1, . . . , 5 corresponds to 3. 



From explicit formula ( 65 ) , it is readily seen that vector "P e q determines a quantum state iff 

i[7>eq(+l/2,ni) +^(-1/2,11!)] + ELi[^eq(+l/2,n fc )-P eq (-l/2,n fc )](l fc • (0,0,1)) > 0, 
\ [Peq(+l/2, ni) + PeqC-1/2, m)] 2 - feLi PV+1/2, n*) - P eq (-l/2, n fc )]lfc) 2 > 0, 



(71) 



i[Peq(+l/2,ni) + 7>cq(-l/2,ni)] - ELi[^eq(+l/2,n fc ) -P eq (-l/2,n fc )](l fe • (0,0,1)) > 0. 



Taking into account the relation 7 3 e q(+l/2, n^) + V e q(— 1/2, n&) = 1/3, k = 1, 2, 3, the obtained system of 
inequalities can be easily depicted (Fig. |4^i). Indeed, using such a constrain on the probabilities, the five- 
simplex for six-component vector "P is identified with the interior of the cube < V e q(+l/2, n^) < 1/3, 
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k = 1, 2, 3. Quantum states are those points on the simplex that satisfy the conditions (71 ). In particular, 
if {n^}? x form an orthonormal basis in M 3 , the quantum states are associated with the ball 

p cq (+l/2,m)-^ + (V eq (+l/2,n 2 )-0 + (V eq (+l/2,n 3 )-^) < ^. (72) 

The general case of arbitrary directions {rife}| =1 is shown in Fig. 4aQ 

In the case of qutrits (j = 1), we restrict ourselves to a numerical solution of the system of seven 



inequalities analogous to (71). It is worth noting that the quantum domain on the fourteen-simplex is 
given by algebraic inequalities — three inequalities of the first order, three inequalities of the second order, 
and one inequality of the third order. Different cut sets of this simplex by hyperplanes V e q(— 1, 114) = 
V eq (—1, 115) = const and P e q(0, nfe) = const, k = 1, . . . , 5 are illustrated in Fig. 4bj^] The cut set of qutrit 
states is a third-degree body of the vector elements V eq (+l,ni), P^+l, 112), and "P eq (+1, 113), with the 
cut set being located between three planes and three second-degree surfaces. 



10 Conclusions 

To conclude, a bijective map of qudit-j states onto single probability vectors "P has been developed. 
In fact, any quantum state is associated with such a probability vector. Quantum states form a convex 
subset on a simplex of possible probability vectors "P, with the boundary of quantum states being the 
(2j + l)-degree body of vector elements V(m, u). Examples of quantum subsets are presented for qubits 
(j = 1/2) and qutrits (j = 1). Components V(rn,u) are fair probabilities, have a clear physical meaning, 
and can be relatively easily measured experimentally. 

To be precise, V{m,Uk) is a joint probability distribution function of two discrete variables - spin 
projection m and unitary rotation Uk from a finite set of rotations {ufc}i."i ■ The number of rotations N u 
is shown to depend on the type of rotations used. Namely, N u = 2 j + 2 if all unitary matrices Uk are 
elements of the group SU(N) with N = 2j + 1, and N u = 4 j + 1 if u k G 577(2) for all k. The latter case 
is considered in detail. The dequantizer operator U^>(m, n^) specifying the direct map p — y Peq and the 
quantizer operator PW(m, n^) specifying the inverse map Poq - > P are presented in the explicit form for 
an arbitrary choice of directions {nfc}^^ 1 . The kernel of the corresponding star-product quantization 
scheme as well as intertwining kernels relating ^-representation and tt>-tomographic representation are 
found. 

A subsidiary problem of the optimum choice of directions {n^}^^ 1 is discussed and partially solved 
for the low spin states, with the optimality implying the minimum relative error ||5/5||2 if errors 57-' in 
the measured probability vector "P are presented. 

Last, but not least, different mappings of density operators onto the probability vectors are unified 
within the concept of the inverse spin-s portrait method. The difference between mappings reduces to a 
particular choice of spin-s portraits implicitly used in these transforms; namely, [12 
portraits, whereas [7j[44] extensively employ spin-j portraits. 

1 Dramatic visualization of qubit states in the probability space is available at Wolfram Demonstrations Project: 
|http: / /demonstrations, wolfram.com /RepresentationOfQubitStatesByProbability Vectors /] 

z The domain of qutrit states in the probability simplex as well as the dynamics of the condition number 
H with regard to an arbitrary choice of directions {rifc}f. =1 is visualized at Wolfram Demonstrations Project: 
http: / /demonstrations. wolfram.com/QutritStatesAsProbability Vectors/ 
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